Diffractive paths for weak localization in quantum billiards 
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We study the weak localization effect in quantum transport through a clean ballistic cavity with 
regular classical dynamics. We address the question which paths account for the suppression of 
conductance through a system where disorder and chaos are absent. By exploiting both quantum 
and semiclassical methods, we unambiguously identify paths that are diffractively backscattered into 
the cavity (when approaching the lead mouths from the cavity interior) to play a key role. Diffractive 
scattering couples transmitted and reflected paths and is thus essential to reproduce the weak- 
localization peak in reflection and the corresponding anti-peak in transmission. A comparison of 
semiclassical calculations featuring these diffractive paths yields good agreement with full quantum 
calculations and experimental data. Our theory provides system-specific predictions for the quantum 
regime of few open lead modes and can be expected to be relevant also for mixed as well as chaotic 
systems. 
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I. INTRODUCTION 

Weak localization (WL), the enhancement of reflec- 
tion by coherent superposition of symmetry-related wave 
components, is an ubiquitous phenomenon of wave trans- 
port through disordered media. Experimental observa- 
tions range from optics ('albedo'), to ultrasound trans- 
mission and seismic wavesi 2 -^^ Very recently, trans- 
port of ultracold atoms through disordered optical lat- 
tices has been proposed as a new candidate for weak 
localization. 5 In mesoscopic physics WL is one of the 
hallmarks for phase-coherent electron transport through 
disordered materials*2±£ Its experimental observation is of 
particular conceptual interest in systems where disorder 
scattering is strongly reduced.—^ In these ballistic sys- 
tems the WL peak in the transport resistance has been 
shown to be sensitive to the underlying classical regular 
or chaotic dynamic s 8 ' 10 : 11 and to the presence of scatter- 
ing resonances^ 

On the most fundamental level, WL probes the particle- 
wave duality in transport. Identification of the relevant 
'paths' in the particle-like dynamics that lead to both 
constructive interference in reflection and to destructive 
interference in transmission continues to pose a challenge. 
In disordered media such a mechanism can be described 
with the help of diagrammatic techniques that allow to 
sum over diffractive paths due to scattering at impu- 
rities or variations of the potential landscaped In the 
ballistic regime of clean quantum dots at low tempera- 
tures where both elastic l e and inelastic mean free paths 
U become large compared to the linear dot dimension, 
l e .i ^> VA (A: area of the dot), the dynamics is co- 
herent and governed by the scattering of the electron at 
the (hard) boundaries of the quantum dot. Thus new 
concepts are required. One of these includes quantum 
corrections to the transport problem in analogy to the 



'diffuson' and 'cooperon' propagators known from the di- 
agrammatic techniques, assuming a static long range po- 
tential within the cavity^ Another concept relates the 
conductance to a sum over classical paths with quantum 
mechanical phases.— In this trajectory-based approach 
partially time-reversed paths of almost equal length but 
a different number of self-intersections were invoked^ 
The are considered to be responsible for WL, and were 
analyzed based on diffractivei 5 - and chaotic^ scattering. 
The latter approach allows to reproduce the WL predic- 
tions of random matrix theory (RMT) and to fulfill the 
current conservation requirement i 14 ' 16 : 17 ' 18 
Identification of the relevant paths for WL in a system- 
specific geometry remains an open problem. Previous 
theories involve an ensemble average over the sample ge- 
ometry, thereby precluding the identification of system- 
specific paths, and an H — > limit, details of which differ 
from each otherj 16 : 17 This complicates a direct compar- 
ison with either experiment or full quantum simulation. 
Moreover, completely chaotic systems are assumed, thus 
preventing applications to quantum dots without disor- 
der and chaotic scattering. For regular billiards clas- 
sical paths alone are insufficient to reproduce quantum 
transport^ To uncover which mechanism for WL is at 
work for this class of systems, we focus in this letter 
on a prototype structure for a ballistic regular cavity: 
the circular quantum billiard with perpendicular leads 
(Fig. [Iji). The circular billiard was shown both exper- 
imentally and numerically to feature a distinctive WL 
effect ££22. It does not contain any of the classical path 
pairs 14 invoked to explain WL in transmission through 
ballistic chaotic cavities. We choose the geometry of the 
billiard and the number of open lead modes N in accord 
with experimental values.— By comparing exact quan- 
tum mechanical calculations^! with semiclassical approx- 
imations that include diffractive paths^ 2 - we are able to 
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FIG. 1: (a) Geometry of circular quantum billiard with radius 
p = \J\ + 4/-7T (in scaled units) and width of the leads d = 
0.25 (~ 0.16p) with one pair of time- reversal symmetric paths 
q and q depicted, (b) Internal diffractive backscattering at the 
lead, described by diffraction amplitude u(#2,#2,fc) (plotted 
is Ke[v(9' 2 ,02 — j , k)] cos(kr) with r as the radial distance 
from the center of the lead and k = ^|2L ) , 



identify the origin of the WL dip in the transmission in 
terms of paths which are diffractively backscattered into 
the cavity and thus are intrinsically cross correlated with 
paths contributing to reflection. 



II. SEMICLASSICAL METHODS AND 
QUANTUM MECHANICAL PATH SPECTRA 

Ballistic transport is described by a wave number (k) 
and magnetic field (B) dependent Hamiltonian S matrix 

with elements Sn,m(k,B) for scattering from mode m in 
lead i to mode n in lead j. For a two-terminal device, the 
transmission amplitudes from lead 1 to lead 2 are denoted 

(2 1) 

as t nm — Sn,m while the reflection amplitudes back into 
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lead 1 are denoted as r. 
given by the Landauer formula 



The conductance is 



N N 



g{k, B) = 2 -^T{k, S) = ^EE I W fc < B )| 2 < C 1 ) 



with T(fc, B) being the total transmission probability. 
The projection of the lead modes 4> m (yi) (yi is the trans- 
verse lead coordinate) onto the constant energy propa- 
gator G(fj,fi,k, B) determines the S matrix elements. 
The transmission amplitudes t nm (k, B), e.g. at B = 0, 
are given by 



tnm{^i B — 0) — i^J k X2 ^ n kxi ,m 

x j dV2 J dy x cj)* n {y 2 )G{ 

r 2j r ii k, B — 0)| a;ii:I ; 2 (f> m 

where k Xltm (k X2tn ) is the longitudinal wave number and 
x\ (x^) is the longitudinal coordinate in lead 1 (lead 
2). Here and in the following, we use atomic units 
(h= \e\ =m eff = 1). 

Standard semiclassical approximations (SCAs) to Eq. [5] 
proceed in two steps. First, the quantum mechanical 
Green's function G(fj,fi,k,B = 0) is replaced by its 



semiclassical limit, 



G^ A (f j ,f f i ,k,B = 0) 
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exp [iS q (k) - i-fiq] , (3) 



which contains a sum over classical paths q(fl — > fj) con- 
necting the leads i and j. The weight of each path is given 
by the deflection factor D q (k) which is a measure for the 
divergence of nearby trajectories. The phase is given by 
the classical action S q (k) which, for a constant potential 
in the interior of the dot, reduces to S q (k) = kL q (L q 
is the length of the path q). In our semiclasssical cal- 
culations, we introduce a small magnetic field B via the 
Aharonov-Bohm phase, S q (k, B) = kL q + Ba q /c (a q the 
directed enclosed area of the path q) . The M aslov index 
[i q is determined by the topology of the classical trajec- 
tory. In the second step, the integrals over the trans- 
verse lead coordinates in Eq. [5] are evaluated approxima- 
tively invoking small de Broglie wavelengths (large k). 
The most common approximation is the stationary phase 
approximation (SPA) which enforces momentum conser- 
vation at the lead cavity junction^ Going beyond the 
SPA, diffraction effects have been introduced on the level 
of the Kirchhoff diffraction approximation, 23 the Fraun- 
hofer diffraction approximation (FDA) ^ or the 'uniform 
theory of diffraction' (UTD) ) 2526 which is based on the 
'geometrical theory of diffraction' (GTD)i22 For all of 
these diffraction approximations the transmission ampli- 
tudes take the following form, 

t S r! A (k,B) = jL=y/k X2in k Xum J2c n (e 2 ,k)J\D q {k)\ 

9 

7T 1 

x exp \iS q (k,B) -i-/igjc m (^i,fc), (4) 

where all injection (ejection) angles 9 E (— 7r/2,7r/2) con- 
tribute. The classical paths q with specular reflections at 
the boundaries get weighted by the diffraction amplitudes 
c m (di,k) [c„(#2, k)] corresponding to the angle Q\ {62) at 
which they enter (exit) the cavity. By improving the ap- 
proximations used for solving the diffraction integrals, an 
improvement of the accuracy of the standard SCA can 
be achieved^ In all semiclassical results presented be- 
low we evaluate the diffraction integrals using the GTD 
(extended by the UTD) to obtain analytical approxi- 
mations for the diffraction amplitudes c m (0,k) featuring 
much better accuracy than within the FDA. In spite of 
this considerable improvement, the standard SCA does 
not reproduce quantum transport properly (even if one 
solves the diffraction integrals numerically exactly, see 
Ref. [l9h . In the following we give qualitative as well as 
quantitative arguments for this failure. As common to all 
standard SCAs, the paths q contributing to the transport 
amplitudes in Eq. ((4]) are purely classical in the interior 
of the cavity connecting the entrance and the exit lead. 
We identify the missing pseudo-paths which are not con- 
tained in Eq. (0| but essential for weak localization. 
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WL appears as a dip in the k averaged transmission 
probability (T)aic and as a peak in the averaged reflec- 
tion probability (R)Ak near B = 0. The full anticor- 
relation, S(R)Ak = —S(T)Ak, highlights the fundamen- 
tal issue of semiclassical descriptions in terms of inter- 
fering paths or, more generally, particle-wave duality in 
WL: Pairs of time reversal symmetric paths, q{r\ — > ri) 
and q{fi — > r*i), contribute to the enhanced reflection, 
8{R)&k , due to constructive interference at B = 0. The 
latter follows from the fact that the phase difference 
S q (k,B) — Sq{k,B) vanishes for time-reversal symme- 
try related pairs as B — ► 0. Such pairs, however, are 
uncorrelated to those pairs of trajectories q'{r\ — ► f^), 
q"(ri — > ^2) which may interfere in transmission. Yet, 
$(T)Ak and S(R)Ak are required to be anticorrelated. 
Trajectories that contribute to transmission should there- 
fore be coupled to those leading to reflection, or in other 
words, subsets of classically disjoint trajectories 'must 
know about each other'. The deterministic outcome of ei- 
ther reflection or transmission must therefore be replaced 
by a probabilistic superposition of both transmission and 
reflection due to particle-wave duality. Correspondingly 
the WL dip in transmission can not be reproduced by em- 
ploying Eq. ([4]) (see Fig. [3J: and discussion below). The 
failure of Eq. (j4|) to properly account for the WL effect 
thus implies that additional quantum effects need to be 
included also for the propagation in the interior of the 
cavity. 

We have identified such non-local quantum correlations 
between trajectories for the circular billiard. Starting 
point is the Fourier transform of the exact quantum am- 
plitudes (for details see Ref[2l|) which include all paths 
of the full Feynman propagator, 

&L,a) = JdkJ dBe^ kL+ ^S^(k,B). (5) 

The Fourier conjugate variables to k and B are the 
length L and the directed area a enclosed by the 
corresponding path. 21 i 22 ' 23 i 24 The resulting length-area 
spectra for both transmission and reflection provide 
unbiased information on the (non) classical properties 
of the whole set of paths of length L and enclosed 
area a entering Feynman's propagator, in particu- 
lar of non-local correlations between reflection and 
transmission without invoking any semiclassical limit. 
The two-dimensional length-area spectra for reflection 
(Fig. [2^ for |f 2 2(i,a)| 2 ) and transmission (Fig. [2b for 
|f 2 2(i,a)| 2 ) of the circular billiard, calculated with the 
help of the modular recursive Green's function method 21 
(MRGM), can directly be compared to paths entering 
the semiclassical approximation. Paths within the 
standard SCA are classical between entering and exiting 
the circular cavity and are marked by dots in Fig. [5] (a) 
and (b). 

Maximal constructive or destructive interference at 
B = requires pairs of paths with (almost) identical 
length (on horizontal lines of fixed L) but symmetrically 
placed on both sides of the a axis. The classical length- 
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FIG. 2: (Color online) 2D path length-area spectrum of the 
exact quantum amplitude \S22{L, a)\ 2 (Eq.[5]) with integration 
ranges k € [2.2n/d, 3.45vr/d] and B/c £ [-3,3] (a) for reflec- 
tion \f22(L,a)\ 2 and (b) for transmission \t22(L,a)\ 2 . Strength 
given by intensity of color, the dots mark classical paths. The 
pairing of classical and peudo-paths contributing to WL is 
visualized [insets in (b)] and examples of diffractive contribu- 
tions by pseudo-paths are highlighted (red markings) . 



area spectrum for reflection (marked by dots in Fig. [2k) 
features such symmetric branches with respect to the 
a = axis and accounts, to a considerable amount, for 
the constructive backscattering. By contrast, classical 
paths contributing to transmission (marked by dots 
in Fig. [5b)) lie on branches on opposite sides of the 
a axis vertically displaced relative to each other and 
thus cannot produce a WL dip. New diffractive paths 
are required to fill the gap and provide approximately 
symmetric branches. Two families of the shortest of 
these diffractive contributions near L ps 10 and 20 are 
highlighted. 

It is now instructive to identify the origin of the paths 
restoring the WL dip in transmission by analyzing the 
pathlength and area distribution. For example, members 
of the shortest non-classical branch are composed of the 
shortest direct classical transmission path qi(f\ —►fa) 
that diffractively backscatters at the exit lead with 
amplitude v(6' 2 ,02,k) (Fig. [lb) followed by a clockwise 
propagating polygon-shaped path (72(^2 - * ^2) that 
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emerges from the exit lead and returns to it [left inset in 
Fig. [2b]- Such paths are referred to as pseudo-paths. 22 
The pseudo-path q(r 1 — > r 2 ) in Fig. [2b (left inset) 
destructively interferes with a classical path q'(r\ — > r*2) 
that first misses the exit lead and subsequently encircles 
the billiard in counterclockwise orientation [right inset 
in Fig. [^h]. (The relative phase between the classical 
and the pseudo-path at B = involves the difference 
in length, the difference in the Maslov indices and the 
phases from diffractive scattering, and is approximately 
7r.) A segment of this second shortest member of the 
family of transmitted pseudo-paths has the shape of 
a triangle [left inset in Fig. [2}d] . Longer paths of the 
same family contain n-polygons which converge in the 
n — > oo limit to creeping trajectories along the circular 
boundary, and determine the endpoint of the branch (see 
Fig. [2]) . This example directly illustrates how diffractive 
wave-scattering at the lead mouth causes the non-local 
coupling of transmitted and reflected trajectories: the 
reflected classical trajectory (72(^2 — > ^2) becomes a 
segment of the pseudo-path q(fi — > r?) contributing to 
transmission. This interplay repeats itself for families 
of longer and more complex paths. The underlying 
mechanism at work here is reminiscent of 'Hikami 
boxes— by which piecewise 'classical' trajectories are 
linked to each other by diffractive scattering events. 
The corresponding pseudo-path semiclassical 
approximation 2 ^ (PSCA) allows to complement the 
present quantum analysis by an explicit summation 
of both the classical paths and the pseudo-paths. 
Within the PSCA diffractive changes of the direction 
in the interior of the cavity are incorporated into the 
Green's propagator in terms of a Dyson equation, 

G PSCA = G SCA + gSCA y gPSCA _ Here ^ unperturbe d 

Green's function G SCA coincides with Eq. [3] including in 
addition paths that are geometrically reflected off the 
leadi2£ The perturbation V represents the diffractive 
scatterings at the lead mouths with scattering vertices 
v(6',6,k) for which we employ the GTD extended by 
UTD. The iterative solution to i th order in V contains 
pseudo-paths consisting of i + 1 segments of classical 
paths interconnected by i diffractive scatterings. The 
corresponding length-area distribution of pseudo-paths 
resulting from PSCA (not shown) closely mirrors the 
respective quantum spectrum thus unambiguously 
establishing pseudo-paths as partners in the path pairing 
required for WL. 



III. RESULTS FOR WEAK LOCALIZATION 

To quantify the role of pseudo-paths in WL we 
have calculated their contribution within PSCA. Since 
only paths up to a finite number of scattering events 
can be numerically summed up and all classical and 
pseudo-paths up to the same length must be included, 
the summation must be truncated for technical reasons 
beyond a maximum pathlength L max . This limitation, 
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FIG. 3: Weak localization (WL) in reflection and transmis- 
sion for the circular quantum dot (see Fig. [TJ with two open 
modes and wave-number average [2.27r/d < k < 2.87r/rf]. (a) 
Quantum calculation of WL with truncated set of quantum 
paths (L < I/max = 40) (see Eq.0. (b) WL calculated within 
the PSCA and the same truncation (no. of paths ~ 7 x 10 6 ). 
(c) WL calculated within SCA (only classical paths) and 
L < L max . Note that a dip in T is completely absent, (d) 
Full quantum calculation (no truncation). 



however, does not impede a quantitative comparison 
with the full quantum calculation since the same 
truncation can be imposed on the true quantum paths 
(see Fig. [2]) by an inversion of the Fourier transform 
(Eq. [5]) with a window for the pathlength < L < L max 
imposed. The resulting PSCA and quantum WL dips in 
transmission (conductance) and WL peaks in reflection 
are in good quantitative agreement with each other 
(Fig. [3k., b). We find satisfactory agreement also on 
the level of the fully differential k-dependence of the 
mode-to-mode scattering amplitudes. 29 Overall, the 
deviation between PSCA and quantum mechanics in 
the A;-averaged mode-to-mode transmission (reflection) 
as a function of the magnetic field is between 0% and 
5% (0% and 10%). Residual differences between the 
PSCA and the quantum results can be attributed to the 
approximations to the amplitudes v(9', 6, k) and c m (9, k) 
involved which can be quantified by direct comparison 
with the exact quantum diffraction amplitudes^ 
To underline the importance of diffractive scattering in 
the interior, we emphasize that the standard semiclas- 
sical approximation without pseudo-paths produces no 
WL dip in transmission (see Fig. [3k). Pseudo-paths (see, 
e.g., the highlighted paths in Fig. [3h.) also considerably 
improve the agreement for the WL peak in reflection 
(compare Fig. [HJd with (3J;). For completeness we also 
show the full quantum calculation for WL without 
truncation (Fig. [3Ji). The influence of longer paths 
beyond the truncation limit (in the present example 
L max = 40 corresponds to 13 traversals through the 
billiard) manifests itself in restoring the unitarity. Note, 
however, that our calculations with truncation reproduce 
the anti-correlation 5(R)Ak — — $(T)Ak, demonstrating 
that anti-correlation between transmission and reflection 
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is an intrinsic property of non-classical path correlations 
and not just a corollary of unitarity. We also emphasize 
that the truncation employed renders the shape of the 
WL peak effectively Lorentzian customarily associ- 
ated with chaotic dynamics. In the limit L max — ► oo 
the triangular shaped peak associated with regular 
dynamic o 10 ' 12 is recovered. This result closely mirrors 
the experiment,' where a transition from a triangular 
to a Lorentzian WL-peak was observed for increasing 
temperatures (there the 'truncation' of paths is imposed 
by the temperature-dependent mean free paths for 
inelastic or dephasing collisions, L ma x = h)- 



IV. CONCLUSIONS 

We have unambiguously identified pseudo-paths which 
are diffractively backscattered from the dot openings to 
be essential for the WL correction in the circular quan- 
tum dot. Our approach neither involves any semiclassical 
limit such as H — > (or N — > oo) nor any assumptions 
about chaoticity or ensemble average. We believe that 
the diffractive contributions identified above provide an 
essential clue to the puzzle on WL at small lead mode 
numbers and in regular systems. The contributions 
result from short paths well below the Ehrcnfcst time 
t/.— — — (unless diffractive corrections arc included 
in the definition of te^)- Our results are in line with 
experiments-^ in which WL is observed for moderate 
values of wavenumbers k with only a few modes TV open 



and a phase-decohering mean free path U extending 
only over a few times the diameter of the structure. 
The latter precludes significant contributions to the WL 
from L > li. We therefore suggest that diffraction plays 
a dominant role for finite h in chaotic systems as well. 
This raises interesting questions regarding the relation 
between our results and previously proposed theories for 
chaotic dots in the H — > limit i 14 i 16 ' 17 i 18 In the present 
case, diffraction originates from the sharp edges at the 
lead. As long as the de Broglie wavelength Xd remains 
larger than the radius r c of the curvature of the rounded 
corners diffraction persists. (In Ref . l3ll it was found that 
the amplitude of WL does not change upon rounding the 
lead edges for small mode numbers.) For Ad <C r c , the 
quasi-classical electron 'sees' a convex-shaped boundary 
and the classically regular billiard turns chaotic. The 
limit h — > thus constitutes the crossover between 
strong contributions from diffraction to emergence of 
chaotic paths scattered at the rounded corners. We 
conclude that the present theory complements rather 
than contradicts previous models for the H — > limit. 
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